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ABSTRACT 

This paper discusses the relationships between the 
coordinates of a point on the surface of an oblate spheroid and 
the coordinates of the projection of that point in several common 
map projections. Since several of the projections are conformal, 
some background material is presented which summarizes the 
theory of conformally mapping ein oblate spheroid to the plane. 
Then, for each projection considered, the equations which map 
the spheroid to the plane and their inverses are gi\ en 
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1. INTRODUCTION AND BASIC DEFINITIONS 

The map projection program MAP2 requires the equations relating 
line and sample in some standard map projection to latitude and longitude 
on the surface of planet. The program must go in both directions; that is, 
given latitude and longitude it must be able to determine line and sample and 
also given line and sample it must be able to determine latitude and longi- 
tude. Equations for line and sample for many different projections from a 
sphere to a plane abound in the literature. The inverse equations are hard 
to find and in some cases are not easy to obtain via algebraic manipulation of 
the direct (line- sample from latitude -longitude) set. Equations in either 
direction for projections from the oblate spheroid to the plane are also rare 
in the literature. 

Reference 1 contains the equations in the inverse direction (obtain 
latitude- longitude from line- sample) as they were used by MAP2 up until 
recently. However, those equations are valid only for projections from a 
sphere, except for the orthographic projection which properly handles the 
spheroidal case. Also Reference 1 contains no direct equations for any 
projection. The purpose of this document is to bring together in one place 
the twelve sets of equations currently used by MAP2. These are the equa- 
tions for the six projections currently implemented; Polar Orthographic, 
Oblique Orthographic, Polar Stereographic, Oblique Stereographic, Two- 
Standard Lambert Conformal Conic, and Mercator projections. Equations 
for both directions are included. The equations are exactly correct for the 
spheroidal case. Only the inverse equations (latitude -longitude from line- 
sample) for the orthographic are identical to those in Reference 1. All the 
other sets are either not in the reference (all direct equations), or have 
been extended to correctly account for oblateness. In addition, for the 
oblique orthographic, the technique given in Reference 1 for the determination 
of which of two solutions to a quadratic equation is the desired one has been 
modified. The old technique failed if the pole appeared in the output 
picture. 

Before the equations are presented, some basic definitions and 
standard notation will be given: 
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Equation of a sphere of radius R: 




I 


( 1 ) 


A standard three>dimensional Cartesian coordinate system is fre- 
quently used throughout this document. Its origin is at the center of the 
planet. Its -tz-axis points to the north pole and its +x-axis points to the 
prime meridian. In such a system the equation of an oblate spheroid of 
equatorial radius and polar radius Rp is: 



( 2 ) 


Throughout this discussion west longitudes will be used. East longi- 
tudes are more natural but for historical reasons MAP2 deals in west 
longitudes. 

Latitude can be defined in several ways and on a sphere these defini- 
tions are all equivalent. On the spheroid, however, they are not. The two 
most commonly used latitudes are geocentric and geodetic. Geocentric is 
defined to be the angle at the center of the planet between the plane of the 
equator and the line to the point of interest. Geodetic latitude is the angle 
between the plane of the equator and the local normal to the surfice of the 
spheroid at the point of interest. For clarity, west longitude will always 
be denoted \ ; geocentric latitude 6; and geodetic latitude 4». MAP2 assumes 
that geocentric latitude is always desired. That is, in the direct mode west 
longitude and geocentric latitude are supplied and line and sample are desired. 
In the inverse mode line and saunple are given and we are asked to compute 
west longitude and geocentric latitude. Howeveri geodetic latitude and another 
type of latitude, called conformal, will frequently be used during intermediate 
calculations. Figure 1 shows the relationship between geocentric and geo- 
detic latitude. Conformal latitude is discussed in Section II. It has no obvious 
geometrical significance on the spheroid. 


2 



77-7 



It i i easy to compute G from 4> or vice versa; 



sign of <)> = sign of 6 . 

If 4> = 0, then 0=0. If <t> = tt/ 2, 0 = w/2. For the remainder of this 
document angles will be specified strictly in radians. 

Two radii are of importance at an arbitrary poinc on a spheroid. One, 
which will be called R, is the distance from the center of the planet to the 
point. 


R(e) 


^eq 




cos^ 0 + sin^ 0 


(4) 


3 



77-7 


The other, which shall be called N, is the local radius of curvature 
of the surface in the direction perpendicular to the local longitude meridiin. 



Note that geocentric latitude enters into the equation for R, whereas geodetic 
latitude is used in the equation for N. The vector from the center of the 
planet to a point at (X, 6 } in our standard 3-D system is: 

X = cos X cos 0 

y = - R (0) sin X cos 0 (6) 

Z = R(0) sin 0 

The minus sign in the equation for Y is due to our use of west 
longitudes. 

Whenever a radius, be it Rp, ^gq» R(®)» N(<J>), appears in an 

equation involving the output (projected) image, it will be assumed to be in 
units of pixels, not kilometers. If F is the scale at the equator (Mercator), 
standard parallels (Lambert), or center of projection (Orthographic and 
Stereographic), in kilometers per pixel, then 

R (pixels) = R (km)/F. (7) 

The 2-dimensional system used to describe the projection plane follows 
MAP2’ s and IPL' s conventions. Its origin is at line 0» sample 0, Its +x-axis points • 
to increasing sample (right); its +z-axis points toward increasing line (down). 

In all projections some longitude will have special significance. In 
Oblique Orthographic and Stereographic it is the longitude of the center of 
projection. In the polar projections it is the longitude which points up 
(decreasing line) as seen at the projection of the pole in the projection plane. 

In the Lambert it is the longitude of the central meridian (which projects to 
a vortical line and is the only such longitude). In the Mercator it is simply 
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.in arbitrary longitude whose sample in the projection is known in advance. 
In all cases the special longitude will be denoted X^. In the oblique pro- 
jections the latitude of the center of projection is 9^ (or 4 >q for geodetic). 
In addition, the oblique projections have associated with them an angle, 
denoted which is the angle in the projection plane of north, measured at 
the center of projection clockwise from up (See Figure Z), 



Figure 2. The North Angle 

In each projection there is some special point whose projected coor- 
dinates are XC, ZC (sample, line respectively). For the oblique projections 
the special point is the center of projection, so (X^, 6 q) projects to (XC, 

ZC). For the Lambert and the polar projections the special point is the one 
(the only one) pole which is visible on the projection. These projections 
have associated with them a variable denoted by CAS. 


CAS = +1 if visible pole is north pole 
CAS = -1 if visible pole is south pole 


( 8 ) 


For the Mercator projection the special point is the point on the 
equator at longitude Xq. Thus the longitude X^ will project to sample XC, 
and the equator will project to line ZC. The various projections are 
characterized by: 

a) Polar Orthographic and Polar Stereographic: F, \q, CAS, 

XC, ZC. 

b) Oblique Orthographic and Stereographic: F, Xq, 6q, ‘I', XC, ZC. 
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c) Lambert: 0^, 0^ (latitudes of standard parallels), F, XC, ZC, 

(CAS is determined by 0^ and 0^ . . . the visible pole has the 
same sign latitude as 0j and 02 ). 

d) Mercator : Kq, F, ZC, XC 

For each projection the characterizing parameters are assumed to be 
known in advance. In MAP2 they are either specified as parameters or, 
if defaulted, are pre-computed by special subroutines which set the param- 
eters in order to center the input picture in the output picture and fill as much 
of the output as possible without loss of data beyond the boundaries of the 
output frame. Where applicable, 4* is determined by minimizing the rotation 
between input and output frames. In fact, it is only for determination of 
default values of the characterizing parameters that MAP2 needs the direct 
equations at all. 

All of the currently implemented projections except Orthographic are 
conformal projections. In order to preserve conformality when the object 
is flattened, a general theory of conformal mapping must be used. Section II 
summarizes the important points of this theory. The results of Section II 
are used repeatedly throughout subsequent sections. Throughout these sections, 
A\ = X. - X^, AX = X - XC, and AZ = Z - ZC. 

II. CONFORMAL MAPPING FROM THE SPHEROID TO THE PLANE 

While it is theoretically possible to conformally project a spheroid 
directly to a plane, the equations are much simpler (especially in the inverse 
direction) if the spheroid is first conformally projected onto a sphere. Then 
the intermediate sphere is projected onto the plane in the usual way. The 
"double" projections so produced differ from the direct projections only to 
terms of fourth order or higher of t , where 
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For Mars this means "accuracy" to one part in 10 if you insist that 
the direct method is better. Since both methods preserve conformality 
exactly, there seems to be no reason to prefer-the direct method, and the 
double projection via an intermediate sphere is used in MAP2, 

In the first projection, from the spheroid to the sphere. Reference 2 
shows that the transformation is given by: 


\' = \ 


<{>' = O' = X , 


( 10 ) 


where the primed coordinates refer to the intermediate sphere (thus the 
equality between geocentric and geodetic latitudes). 

Thus the first projection preserves longitude, and latitude on the 
sphere equals X, the so-called conformal latitude on the spheroid (this is the 
third type of latitude mentioned in the previous section). Reference 2 gives 
the relationship between conformal and geodetic latitudes as: 


where £ is given by Eq, 9. Throughout the rest of this document the Greek 
letter CHI (X) is consistently used for conformal latitude. Obviously, 
if is known, X is easy to compute. If the right side of Eq. 11 is denoted 
by Q, then: 


X 




arctan Q - 



( 12 ) 


sin X 


cos X 


- 1 
+ 1 


2Q 


+ 1 


(13) 


(14) 
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Unfortunately the inverse problem, obtaining given X, is not so 
simple. The equation is transcendental and cannot be solved explicitly for 
4>. MAP2 solves the problem numerically via a successive approximation 
technique as follows: 

First rearrange Eq. 11: 


tan 


(M) = 


72 


(15) 


As a first guess, replace all appearances of 4> on the right with X. Thus 


^/2 


<(> j = 2^arctan Qj - ^ ) 

Now use ^ ^ as a guess for <t> and compute 

«.■ »(M) (;;4^)" 


(16) 


(17) 


18) 


4*2 = 2^arctan Q 2 ‘ T ) * 


( 19 ) 


Repeat this process until 14*., . - 4* • I i some tolerance value, which 

mj ■ 1 X I 

MAP2 has currently set at 10 radians. The process converges extremely 
rapidly - on Mars, where £ * 0. 1, 4 iterations are almost always sufficient. 

An alternative which will converge even faster in some cases is to use 
a Newton’s method technique. This seems especially attractive since the 
derivative of the right side of Equation 11 has a surprisingly simple form, 
namely 


d 


[-( 14 )] 


(1 - £^) sec 4* 
1 - £ Sin 9 


( 20 ) 
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However, sed practical experience shows that this method fails near 
the pole due to the fact that the second derivative of tan (ir/4 + X/2) is 
approximately sec ^ times the first derivative. Near ir/2, sec ^ approaches 
infinity so the basic assumption of Newton’s method, that the function is 
approximately linear, fails. The successive approximation technique of 
Eqs. 15- 19 seems stable at all latitudes. 

The generally Inferior Newton's method technique is mentioned here 
only because MAP2 actually uses it in the Lambert conformal conic pro- 
jection. The code for this projection was the first to be written, before the 
method's drawbacks were discovered. Since this projection is rarely 
applied near the pole, it was not felt desirable to modify an already working 
code. 

To summarize, all conformal projections from a spheroid are easily 
handled by using the standard spherical equations with two easily imple- 
mented modifications: 

1) Replace all latitudes in the spherical equations with X, the 
conformal latitude on the spheroid. 

2) For the radius of the sphere in the equations use the radius 

of the intermediate sphere. This radius, is given in 

Reference 2 as 


R. . 
int 


'I’o 

CCS Xrt 


( 21 ) 


where the zero subscript refers to the center of projection. 
Equation 21 breaks down at 4> = 0 = X = 0 and at <|> = 0 = X = ir/2. So 


-f 



( 22 ) 


The remainder of this document is a description of the equations on a 
projection by projection basis. For all except the orthographic, the 
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techniques of this section make the problem conceptually simple, though 
algebraically complex. For the orthographic projection, which is not con- 
formal, Section II does not apply. But the simple geometrical definition of 
the orthographic projection saves us and allows the problem to be solved. 


UI. POLAR ORTHOGRAPHIC 


A, Direct Equations 

On the spheroid we define the orthographic projection to be one of true 
perspective. The plane of the projection is tangent to the spheroid at 
the center of projection. The perspective point is at infinity, so the 
perspective lines are all parallel to each other, perpendicular to the plane 
of projectit’i, and thus parallel to the local surface normal at the center of 
projection. 

The oblateness of the spheroid hardly raises any complications in the 
polar case, although it causes considerable trouble in the oblique case. In 
the polar case longitude meridians are straight lines radiating from the pole 
at their true orientation. Latitude lines are circles having their true 
radius. Thus we have 


X = R(0) cos 0 sin ’S' CAS + XC 
Z = - R(0) cos 0 cos AX. + ZC 


(23) 


B. Inverse Equations 

The above equations are easily inverted algebraically to produce 


RAD = 



+ AZ' 


(24) 


A.V 


arccos 


AZ \ 
RAD/ 


(25) 
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X = CAS * AX - Aq 


r| Rp 

RAD^ 

eq 

If the qusintity under the radical is negative, the point specified is not 
on the planet. Remember that the orthographic projection does not fill the 
projection plane. 


9 = CAS arctanl 



IV. OBLIQUE ORTHOGRAPHIC 


A. Direct Equations 

Consider Figure 3. 



Figure 3. Oblique Orthographic Projection on the Spheroid 

First compute the vector from the planet center to the point of interest in 
our standard 3-D coordinate system. 

X = R(9) cos X cos 6 

Y = - R(0) sin X cos 6 (27) 

Z = R(0) sin 0 
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Next rotate these coordinates to a system whose z-axls is parallel 
to the perspective lines, and whose fx-axis points south (in the projection 
plane, which is now parallel to the x-y plane). It should be clear that in 
the new coordinate system the x and y coordinates are directly related to 
line and sample in the projection plane. The required rotation is a double 
one. First rotate by - about the z-axis; then by tt/2 - 4*q about the new 
y-axis. Note that the geodetic latitude of the center of projection is involved 
here. Thus we obtain 




\ 



/ sin <|>Q cos Xp 


sin X 


^cos 4 >q 


0 


cos X 


0 


- sin <j> Q sin Xq 

- cos <j>Q \ . 

lA 

cos Xq 

0 

Y 

- cos <t> Q sin Xq 

»in+o / 

/ 


(28) 


Multiplying out Eq. 28, substituting Eq. 27 into the result, ignoring Z', 
algebraically simplifying, and renaming variables to conform more closely 
to the MAP2 2-D standard X - Z system, results in 


X" = - R(0) cos 6 sin AX 

Z" = R(0) [sin C03 0 cos AX - cos sin 0] 


(29) 


The double primes in (29) are there as reminders that we have not yet 
reached our standard system in which X points to increasing sample, Z to 
increasing line, with the center at (0, 0). Due to the rotation in (28), Z" 
points south, not to increasing line. Also the origin is not (0, 0) or even the 
center of projection, but instead is the intersection of the projection plane 
and the perspective line passing through the center of the planet (point Q 
in Fig. 3). First, translate so that the center of projection is at the origin: 


X' = X" 


Z' = Z" - R(0q) sin (<|)q - Oq) 


(30) 
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Finally, rotate about an axis perpendicular to the projection plane by 
-4' and translate to (0, 0) as origin: 


X = X' cos - Z' sin 4> + XC 
Z = X' sin 4^ + Z' cos 4^ + ZC 


(31) 


Equations 29, 30, and 31 suffice to compute X and Z for any specified 
(\, d) in the oblique orthographic projection. 


B. Inverse Equations 

These equations are presented here in cookbook style for the sake of 
completeness. They were originally derived by Avnie Schwarz and are 
described in detail in Reference 1. The basic plan is to determine, using 
analytic geometry, the point of intersection between the spheroid and the 
perspective line passing through the given point on the plane of projection. 


n = 0Q + arccos 


cos^ 0. 


SR. 


sin^ 0. 


-- + 


~Z ! 2 

cos 0Q sin 0Q 


R 


eq 


(32) 


L = 



(33) 


6 = 4^ + arctan 


-AX 

-AZ 


(Remember, 4< = north angle) 


(34) 


DD = J (AX)^ + (AZ;^ 


(35) 
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A1 

B1 

Cl 

D1 


- sin 6 sin Xq - cos 6 cos Xq sin 

- sin 6 cos Xq + cos 6 sin Xq sin ^ 

cos 6 cos n 
2 

-DD sin 6 + R(®q) cos 6 sin U cos 9 q - ^0 

2 2 2 2 
-DD cos 6 cos Xq sin 9 q sin fz - DD cos 6 cos Q. 

2 2 2 
-DD sin 6 cos 6 sin Xq 


cos cos & 
(36) 


A2 = - cos 6 sin Xq cos L + sin 6 sin 9q cos Xq 
B2 = - cos 6 cos Xq cos L - sin 6 sin 0 q sin Xq 
C2 = - cos 0Q sin 6 


A2 • 

Cl 

- A1 • 

C2 

A2 • 

B1 

- A1 • 

B2 

B1 • 

C2 

- B2 • 

Cl 

- Y 





(37) 


G2 = 

eq ' 

DRPSQ = Rp 


(38) 


Z1 = DRPSQ + Y^) + G2 • p ^ 


,,, o- C2 • D1 • DRPSQ + 3 • B2 • D1 • G2 
Ki = r 


= Y\ 


(39) 
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K3 = [zo • C2 • p • B2 • - (C2)^ DRPSQ - (B2)^G2 - a^(B2)^R^^ 

- P^(C2)^R^ ](D1)^ + G2* DRPSQ+R^ (DRPSQ)©^ 

/ 

+ R^ p^ • G2 
eq ^ 


At this point check to see if K3 <0. If it is, the point originally 
specified is off the planet - there is no intersection between the spheroid 
and the perspective line. 


XI = K1 + K2 -/kT 
X2 = K1 - K2 ViO“ 

Yj _ -D1 • C2 + tt- XI 


Y2 = 


Z1 = 


-D1 • 

Y 

C2 + o • 

X2 


Y 


-B2 • 

D1 + P • 

XI 


e 


-B2 • 

D1 + P ' 

X2 


(40) 


You now have two vectors, (XI, Yl, Zl) and X2, Y2, Z2), which are 
the lines from the planet center to the two points of intersection between 
the spheroid and the perspective line, expressed in our standard 3-D 
coordinate system. Call the vectors Vj and respectively. 

Next compute the vector P from the planet center to the center of 
projection 


P^ = R(0q) cos 0q cos \q 
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?Y = - R (6q) cos sin (41) 

Pz = ®0 



Choose the vector from Vj and which minimizes the quantity 
- P I and - pj. Call the winning vector's components X, Y, and Z. 


X = 


- arctan 



(42) 


0 = 


arctan 



V. POLAR STEREOGRAPHIC 


A. Direct Equations 

The equations in Subsection A are taken directly from Reference 3, 
which used the principles of Section II. First compute p, the length in 
pixels on the projection of the line from the pole to the point of interest. 


p = 2 R (1 + e (1 _ £ j-(l+p-)/2 

eq 


_ CAS»<(> I 


e/2 

/ 1 + e sin (CAS*4>) \ 

\T'^~e"TIn'(UAS^^^ ) 


(43) 


X = CAS • p sin AX + XC 
Z = - p* cosAX + ZC 


(44) 
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B. Inverse Equations 

These equations follow from simple algebraic manipulation of the 
equations of Subsection A, above. 

First, compute 

C= 2R ( 1 + (1 - 

eq 


and 


P = 


J 


(AX)^ + (AZ)^ 


(46) 


Then 


X = arctan ( ' ) + X^ (47) 

cIs [7-»"‘“(fe)] (■‘81 

Use the equations of Sections I and II to obtain 4> and 6 from X. 


VI. OBLIQUE STEREOGRAPHIC 


A. Direct Equations 

I 

The spherical equations were taken from Reference 3 and modified 
according to the principles of Section II, First compute line and sample in a 
coordinate system centered at the center of projection and which has north 
at the top. 


X' 


2R. . cos X sin AX 
int 


1 + sin Xq sin X + cos X^ cos X cos AX 


(49) 
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- 2 £cob Xq sin X - sin X^ cos X cos A\] 

^ “ 1 + sin Xq sin X 4 cos X^ cos X cos 

Now rotate so north is at the desired angle and translate to the usual 
origin: 


X = X» cos 4* - Z» sin 4^ + XC (SO) 

Z = X' sin 4^ + Z' cos 4^ + ZC 


B. Inverse Equations 

Direct algebraic inversion of the equations of Subsection A is very dif- 
ficult so explicit advantage of the existence of the intermediate sphere is taken. 
Analytic geometry is used to find the intersection of the given perspective line 
of interest and the intermediate sphere. The latitude of the point of inter- 
section on the sphere is the conformal latitude of the point of interest on the 
spheroid. The method works because the stereographic projection from a 
sphere is a true perspective projection with plane of projection tangent to the 
Sphere and perspective point on the opposite side of the sphere from the center 
of projection, at coordinates (\q + Xq). See Figure 4. 



Figure 4. Oblique Stereographic Projection of a Sphere 
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Phase I — Translate the given projection plane coordinates to a new 
origin at (XC, ZC) and rotate so north is at the top (X' points east, Z' south) 
from the origin. For ease in conversion to our standard 3-D system, the Z' 
axis will be renamed Y*. This rotation is by 


X' = AX cos 4< + AZ sin ‘I' 
Y' = - AX sin 4* + cos 4^ 


(51) 


Phase 2 — Now consider X', Y’ to be two of the three coordinates in a 
non-standard 3-dimensional system whose origin is the center of the inter- 
mediate sphere. From the definition of X' and Y' above, and in order to 
obtain a right-handed system, the Z* axis of this system points to the 
anticenter of projection (see Figure 4). 

Thus, the coordinates of the projection of the point of interest in this 
new system are: 


Xj = X' 

Yj = Y' 

Z, = - R. ^ 
1 mt 


(52) 


Similarly the coordinates of the anti-C. P. are; 


= 0 
= 0 


(53) 


We wish to find the intersection of the perspective line of interest, 
which passes through (Xj, Yj, Z^) and (X^, Y^, Z^), and the intermediate 
sphere. Equation of the sphere is 




( 54 ) 
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The equation of a line in 3- space is 


X - Xj ‘ '^1 Z - Z 



or, in our case, 


X - X' Y - Y* ^Int 

' -Y' " 2 

int 


( 55 ) 


(56) 


There are two intersections of the perspective line and the sphere. 
One is just the anti-C. P. and is not of interest. Simple but tedious algebra 
shows that the coordinates of the interesting intersection are: 


Z = R. 


xnt 


m2 


+ Y' 


X»2 + V'^ + 4 R 


2 1 
int J 
F^“ 

int_ 


X 


2 R 


int 


(Z . R.^^) 4 X- 


(57) 


Y 


-Y' 


2 R, 


int 


(Z 4 R.^^) 4 Y- 


Phase 3 — Rotate the coordinates just found to our standard 3-D sys- 
tem with 4z-axis pointing to the north pole. This involves first a rotation 
about the x-axis by xr/2 + then about the new z-axis by \q - Thus 

Xp = X sin 4 Y sin Xq cos - Z cos Xq cos X^ 


Yp = X cos Xq - Y sin sin X^ -r Z cos Xq sin X^ (58) 

Z_ = - Y cos X^ - Z sin X^ 

F 0 0 
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(Xp» Yp, Zp) are o\ir desired coordinates, and the latitude and 
longitude follow immediately. 




X = 



( 59 ) 


As usual, ^ and 6 must then be computed from X. To summarize, the 
given X and Z plus the characterizing parameters of the projection allow 
Eq. 51 to be used to determine X', Y'. Then Eq. 57 determines X, Y, Z. 
Equation 58 determines Xp, Yp, Zp and finally Eq. 59 gives X and X. 


VII. TWO-STANDARD LAMBERT CONFORMAL CONIC PROJECTION 


A. Direct Equations 

These equations come directly from Ref. 3. lirst compute K, the 
constant of the cone, as follows: 


ln[N^ »ln^N^eo,»^j 


f (i - ‘ •>" (CAS« ♦ ^ 


- u 

P . t lin (CAS» \ 


(_ \ 1 • ■ lin (Ca5* ♦ ] ^ J 

\ i /. 

\ 1 - 1 lin (CAS' ♦ / 

1 2 


(60) 


where subscripts 1 and 2 refer to the latitudes of the standard parallels. 

Next, compute C, the distance in pixels on the projection from the 
visible pole' to any point on the equator. 


C = 


N j cos ^ ^ 


K 


' /l + t sin (CAS« j) Jl CAS* 4> j \ ' 

_\1 - t Bin (CAS“* ^ j) / \4 2 / ^ 


K 


(61) 


21 


77-7 


Next, get p = d' 'ance on projection from visible pole to point of 
interest: 


P ^ 


/ 1 + E sin (CAS«4>) 

\ 1 - e sin (CA^«<|>) j 


tan 




Finally: 


( 62 ) 


X = - p sin (K • A\) + XC 
Z = CAS* p cos (K • A\) + ZC 


(63) 


B. Inverse Equations 

These follow from direct algebraic manipulation of the equations of 
Subsection A. 

p = J(AX)^ + (AZ)^ . (64) 

If p = 0, you are at the visible pole. 

Check to see if 


Z 


arctan 


AX \ 
CAS*AZ / 
K 


(65) 


has an absolute value greater than x. 

If it does, the point of interest is not on the planet (the Lambert pro 
jection does not fill the plane). If it does not, then 


\ = 



cot 



CAS*x j 



( 66 ) 


Compute 0 from '' in the usual way. 
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/ 

/ 


VIII. THE MERCATOR PROJECTION 


A. Direct Equations 

These were obtained from Reference 3. 


X=-R *AK+XC 

eq 


r 

e/2 


I-) 

• 

Z = - R In 


f T + 

+ ZC 

eq 1 


2/. 

. 


(67) 


B, Inverse Equations 

These are obtained by inversion of Eq. 67: 


\ = 


AX 


R 


+ K 


eq 


( 68 ) 


y ^ V -AZ/R 
tan(M)=e 

Compute 0 from X in the usual way. 

Care must be taken in the Mercator projection to insure that the pro- 
jection, which does not fill the plane, is not accidentally repeated indefinitely 
by ignoring this fact. The absolute value of AX must never be allowed to 
exceed ir. Thus in Eq, 67, if |^" ^q| = ! is greater than x, add or sub- 

tract Ztt to it to ensure that |aX| £ it. In Eq. 6 8, if the quantity )AX |/R 

eq 

> IT, the point specified is off the planet. If these precautions are taken, 

the longitude becomes the central meridian, with exactly half the planet 

on its left and half on its right. Alternatively, you could restrict AX to 

the range - Ztt £ AX - 0 in Eq. 67, and have the off-planet criterion for 

Eq. 68 be: XC = 1 at all times, AX < 0 or AX > Zir R , If this route is 
^ eq 
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taken, \q becomes the longi»:ude at the extreme left edge of the picture 
(sample 1). This latter technique is used by MAP2. 
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